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LINEAR  THEORY  OF  A  LARGE  VOLUME,  HIGH  POWER  GYRO 
TRAVELING  WAVE  AMPLIFIER 


1 .  Introduction 

This  paper  discusses  a  new  concept  for  a  high  power  millimeter 
wave  gyro  amplifier.  While  it  is  potentially  viable  at  any  short 
wavelength,  we  particularly  focus  on  94  GHz,  the  clearest  part  of  the 
atmospheric  propagation  window.  Also  we  envision  this  concept  as 
a  rf  source  for  a  large,  fixed  facility  such  as  the  Haystack  or 
Kwajalein  Radar  facility,  each  managed  by  Lincoln  Laboratory.  Both 
of  these  antennas  are  qualified  for  radiation  at  94  GHz.  The  former 
could  accept  an  rf  source  of  the  size  we  propose  now.  For  the  latter, 
some  reengineering  of  the  facility  would  have  to  be  done} 

The  basic  source  concept  is  shown  in  Fig.(l),  and  a  modified 
version  is  shown  in  Fig.(2).  Instead  of  a  waveguide  configuration,  it 
utilizes  Gaussian  radiation  beams  directly.  This  is  advantageous, 
particularly  for  the  projects  of  interest  to  Lincoln  Laboratory.  They 
have  developed  a  number  of  high  power  components  for  Gaussian 
beams  which  we  will  discuss  shortly.  While  their  inventions  were 
intended  for  the  beam  transmission  system,  there  is  no  reason  they 
cannot  be  used  directly  in  the  amplifier  itself. 

As  an  alternative  configuration  one  could  also  use  a  corrugated 
waveguide  with  mitre  bend  reflectors,  shown  in  Fig.  3,  or  with  U 
bends  (not  shown).  This  corrugated  waveguide  supports  TEM  modes 
having  very  low  loss.  The  waveguide  is  shown  with  large  gaps  near 
the  point  of  electron  beam  traversal.  These  large  gaps  are  to 
eliminate  cutoff  modes  of  the  waveguides  at  these  locations,  as  these 
would  be  the  most  likely  competing  modes.  Since  the  transverse 
mode  structure  in  both  the  quasi-optical  and  corrugated  waveguide 
configurations  are  each  nearly  Gaussian  in  structure,  we  analyze  here 
the  beam  interaction  with  Gaussian  modes. 

The  configuration  in  Fig.(l)  has  a  radiation  field  bouncing 
between  mirrors  in  a  serpentine  path  going  from  right  to  left.  A 
gyrotron  electron  beam  (i.e.  one  with  transverse  energy)  propagates 
through  the  center  of  the  radiation  spot  on  the  vertical  radiation  legs 
in  the  figure.  As  it  propagates,  it  gives  energy  to  the  fields  in  a 
modified  version  of  the  gyro  traveling  wave  amplifier  (TWA) 
interaction.  This  configuration  was  motivated  by  earlier  work  at  NRL 
on  a  gyroamplifier  in  a  folded  waveguide  system^,  shown 
schematically  in  Fig.  (4).  However  in  the  waveguide  system,  even  at 
35  GHz,  the  waveguide  is  quite  small  in  the  transverse  direction. 


Manuscript  approved  March  18,  1996. 


1 


thereby  limiting  the  power  capability.  Our  configuration  is  also 
motivated  by  other  work  at  NRL  on  the  quasi-optical  gyrotron^.^, 
shown  schematically  in  Fig.(5).  Here  the  gyrotron  cavity  was  a  Fabry 
Periot  resonator  oriented  perpendicular  to  the  magnetic  field.  As  the 
beam  traverses  the  resonator,  iv  gives  energy  to  the  fields  via  the 
conventional  gyrotron  interaction.  Recently  it  was  also  pointed  out 
that  the  MIT  group  has  also  analyzed  a  configuration  similar  to  this. 
Their  analysis  relied  mostly  on  particle  simulation.  This  work  here 
complements  that  by  offering  analytic  theory. 

Experience  in  running  this  device  over  about  5  years  has 
shown  that  it  operates  basically  as  a  gyrotron,  but  at  somewhat  less 
efficiency6.7.  This  can  be  partially  compensated  by  the  relative  ease 
of  incorporating  a  depressed  collector,  due  to  the  fact  that  the  beam 
and  radiation  propagate  orthogonally.  The  best  results  of  the 
experiments  so  far  show  an  efficiency  of  20%,  which  has  been 
increased  to  30%  with  a  very  simple  single  stage  depressed  collector. 
The  quasi-optical  gyrotron  has  a  high  axial  mode  density,  but  only  a 
single  transverse  mode  has  ever  been  observed.  Recent  work  has 
shown  a  number  of  ways  in  which  the  device  can  operate  efficiently 
in  a  preselected  axial  mode.  One  of  the  ways  has  been  to  use  a 
prebunching  cavity.  This  prebuncing  cavity  is  also  shown  in  Fig.(5). 

It  was  powered  by  an  extended  interaction  oscillator  (EIO)  at  85  GHz. 
The  prebunched  beam  either  excited  the  main  cavity  as  an  amplifier 
for  the  case  where  the  main  cavity  was  stable,  or  else  controlled  it  as 
a  phase  locked  oscillator  for  the  case  where  the  main  cavity  was 
unstable,  but  not  very  far  above  threshold.  This  amplifier  should  be 
capable  of  high  average  power,  since  the  size  of  all  components  are 

large. 

It  is  important  that  amplifier  operation  has  been  demonstrated 
on  the  quasi-optical  gyroklystron.  Nevertheless  the  device  is  not 
suitable  as  an  amplifier  if  reasonable  bandwidth  is  desired.  Since  the 
cavity  mode  has  very  high  Q  the  bandwidth  was  only  a  few  MHz. 
However  with  the  experience'  obtained  in  the  quasi-optical  gyrotron 
and  gyroklystron,  it  is  natural  to  propose  quasi-optical  gyrotron 
traveling  wave  amplifier  for  the  case  where  reasonable  bandwidth  is 

desired. 

The  configurations  shown  in  Figs.  (1  and  2)  appear  to  have 
several  advantages  over  waveguide  configurations.  First  of  all,  the 
mode  density  is  low.  Assuming  that  the  region  outside  of  the  ray 
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path  uses  millimeter  wave  absorbers,  the  only  allowed  modes  are  the 
modes  of  interest  themselves.  We  envision  a  forward  wave 
interaction.  Since  this  mode  cannot  be  absolutely  unstable,  it  cannot 
oscillate  without  a  source  to  drive  it.  However  the  backwards  mode 
can  also  propagate,  and  this  might  be  absolutely  unstable  depending 
on  the  current  and  length  of  the  system.  Typically  in  such  a  case, 
one  either  designs  the  amplifier  short  enough  to  be  below  oscillation 
threshold,  or  else  uses  severs  in  the  circuit.  As  we  will  see,  another 
possibility  in  this  configuration  is  the  use  of  a  high  power  Gaussian 
mode  circulator,  developed  at  Lincoln  Laboratory. 

A  particular  advantage  of  the  ccmfiguration  proposed  is  that 
there  is  no  zero  group  velocity  mode.  In  traveling  wave  amplifiers,  it 
is  almost  always  the  cutoff  mode  that  interferes  with  operation  by 
breaking  into  spontaneous  oscillation.  In  the  corrugated  waveguide 
configuration  of  Fig.  (3),  the  portion  of  the  guide  in  which  a  cutoff 
mode  is  most  likely  to  be  excited,  has  been  removed.  Furthermore, 
we  will  also  see  that  for  electron  beam  energies  of  interest,  if  the 
amplifier  is  stabilized  at  a  cyclotron  harmonic  n,  it  is  also  stabilized 
for  all  other  harmonics  p  as  long  as  p>n.  Finally,  it  is  clear  that  in  the 
quasi-optical  gyro  TWA,  the  sizes  of  all  components  are  relatively 
large,  meaning  that  the  power  handling  capability  is  also.  The 
development  of  such  an  amplifier  with  peak  power  of  100  kW  and 
average  power  in  excess  of  10  kW  appears  to  be  achievable. 

Let  us  now  continue  by  discusing  two  inventions  developed  at 
Lincoln  Laboratory8.9  which  could  render  simpler  the  development 
of  the  quasi-optical  gyro  TWA.  The  first  is  the  clam  shell  reflector 
shown  in  Fig.  (6).  When  a  Gaussian  beam  reflects  obliquely  from  a 
spherical  mirror,  there  is  some  depolarization.  By  utilizing  the 
double  mirror  system  shown  (the  clam  shell),  the  depolarization  in 
the  first  mirror  is  reversed  in  the  second,  so  the  polarization  is 
preserved.  Each  mirror  in  Fig.  (1)  might  be  a  clamshell,  alternately, 
if  the  clam  shell  mirrors  can  be  separated,  the  two  adjacent  mirrors 
in  each  horizontal  leg  can  be  a  clam  half  of  a  clam  shell. 

The  second  is  the  quasi-optical  millimeter  wave  circulator 
shown  in  Fig.  (7).  It  utilizes  a  ferrite  to  rotate  the  polarization  of  the 
radiation  45°.  Typically  a  permanent  magnet  is  used  with  a  field  just 
large  enough  to  drive  the  ferrite  into  saturation.  Here  the  losses  in 
the  ferrite  are  minimized.  However  it  is  also  possible  to  use  the  main 
field  of  the  gyro  TWA,  although  this  is  a  much  larger  field  than 
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would  ordinarily  be  used.  These  ferrites  are  typically  very  lossy  at 
the  cyclotron  resonance,  but  not  so  lossy  at  other  frequencies  . 
Particularly,  if  the  radiation  is  at  twice  the  cyclotron  frequency,  the 
losses  in  the  ferrite  are  not  that  large.  The  circulator  can  then  be 
used  in  the  circuit,  for  instance  along  the  horizontal  legs  of  Fig.  (1)  to 
prevent  the  backward  wave  from  propagating  through  the  system. 

It  can  be  used  as  an  alternative  to  the  sever,  but  with  the  advantage 
that  one  does  not  have  to  worry  about  dumping  the  partially 
amplified  radiation  somewhere  in  the  tube  (ie  at  the  sever). 
Furthermore,  if  the  gyro  TWA  operates  at  the  harmonic,  the 
circulator  can  not  only  decouple  the  backward  wave,  but  also  provide 
strong  dissipation  at  the  fundamentak _ 

This  paper  analyzes  the  linear  theory  of  the  quasi-optical  gyro 
TWA.  In  Section  2,  the  general  formalism  is  introduced.  Section  3 
calculates  what  is  essentially  E*  J  for  the  configuration.  This  quantity 
is  needed  to  calculate  the  linear  dispersion  relation.  Section  4 
calculates  the  instability  condition.  Basically  this  calculates  the 
maximum  length,  either  for  the  entire  system  or  between  severs  or 
circulators.  Section  5  gives  a  few  examples  and  results.  In  a  sense 
the  quasi-optical  gyro  TWA  is  the  complement  of  the  waveguide 
system.  Where  the  latter  is  fighting  large  sizes,  the  former  is  fighting 
large  size.  Typically  we  find  system  lengths  of  half  a  meter  to  a 
meter  and  magnet  bores  of  4  to  6  inches.  The  appendix  shows  that 
in  the  linear  regime,  the  klystron  effect,  that  is  the  bunching  of  the 
beam  in  one  vertical  leg,  followed  the  extraction  in  the  next,  is 
negligibly  small. 
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2.  General  Formulation 

In  Fig.  (1),  let  the  plane  of  the  paper  be  the  yz  plane  with  the 
constant  magnetic  field  being  in  the  z  direction.  The  natural 
polarization  for  the  electromagnetic  wave  in  a  gyrotron  configuration 
is  then  the  x  direction.  Maxwell’s  equations  for  this  component 
reduce  to 

VxVxExix  +(l/c)^ix3^Ex/3t2  =  -ix(47t/c2)9Jx/dt  (1) 

The  configuration  as  shown  in  Figs.(l)  is  periodic  with  periodicity 
length  L  in  the  z  direction,  and  confined  radially  to  be  within  the 
outer  conductor.  Therefore, by  Floquet’s  theorem,  all  quantities  (E,J 

etc)  can  be  expressed  as  a  periodic  function  of  z  times  an  exponential 
exp  ikz.  Further,  we  imagine  that  the  ray  bounces  around  in  the 
periodicity  cell  along  some  ray  path,  the  distance  along  which  is 
denoted  by  1.  If  we  assume  the  wave  propagates  strictly  in  one 
direction,  say  to  the  right  in  Fig.  (1),  neglect  diffraction,  and  the 
coupling  to  the  backwards  wave,  Eq.  (1)  may  be  reduced  to  a  single 
first  order  scalar  equation 

aEx/at  +  caEx/ai  =  -27iJx  (2) 


and 


Ex=  P(y,z)  exp  i(kz-a)t)  (3) 

where  P(y,z)  is  periodic  in  z  with  periodicity  length  L.  Now  say  that 
the  actual  ray  path  the  radiation  traverses  within  a  periodicity  cell 
has  length  A.  Let  rj.  be  the  spatial  variable  describing  the  radiation 
profile  perpendicular  to  1.  Then  let  us  take  as  an  approximation  for 
P,  the  function 

P  =  R(ri)  expi[kLl/A  -  kz]  (4) 

Note  that  as  the  wave  propagates  a  distance  A,  z  advances  a  distance 
L  so  that  P  is  periodic  in  z  with  periodicity  length  L.  For  the  case  we 
will  consider  here,  R  is  a  Gaussian  with  spot  size  ro, 

R(ri)  =  exp  -(r/ro)2  (5) 
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Let  us  now  stipulate  that  the  electric  field  Ex  can  be 
decomposed  into  a  summation  of  orthogonal  functions 

Ex=SiaiPi  (6) 

For  the  normalization  we  take 

Jd3r  Pi*Pj  =  WA  (7) 

where  in  the  three  dimensional  spatial  integral  above,  the  z  integral 
is  taken  only  over  a  length  L.  Let  us  further  take  Pi  as  being  given 
by  Eq.(4).  If  we  assume  that  in  the  decomposition  of  the  electric 
field,  ai  is  nearly  unity  and  all  other  ai’s  are  small,  we  can  multiply 

Eq.(2)  by  Pi*,  integrate  over  space  and  find 

L 

-i(co  -ckL/A)  =  (itro2A)-l27tlodzJd2r  JxPi*  (7) 

This  then  constitutes  the  dispersion  relation.  To  complete  the 
calculation,  we  must  compute  the  right  hand  side,  which  is  done  in 
the  next  section.  Notice  that  in  the  limit  of  zero  current,  the  local 
wave  number,  kL/A  is  smaller  than  the  wave  number  k 
characterizing  the  Floquet  solution.  This  is  so  that  when  z  chanps  by 
L,  or  1  changes  by  A  (which  is  longer  than  L),  the  phase  change  is  the 
same.  Also  note  that  because  the  k  value  is  much  larger  than  2;t/L, 
it  is  obviously  not  restricted  to  the  first  Brillouin  zone  (indeed  there 
is  no  reason  for  the  k  to  be  in  any  particular  zone).  The  index  of  the 
Brilloun  zone  is  given  approximately  by  (oA/27ic. 
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3.  Calculation  of  /Pi*Jx 


To  proceed,  we  calculate  JPi*Jxfor  the  configuration.  We  will 
calculate  this  only  for  the  configuration  of  Fig.(l),  and  later  will 
simply  state  the  results  for  Fig.(2).  In  each  periodicity  length  the 
electron  beam  the  intersects  the  radiation  twice,  once  going  up,  and 
once  going  dowft.  In  each  vertical  leg,  the  electric  field  is  given  by 

E  =  EoixCxp-(r/ro)^expi((-l)”*ky-a)t+<hin)  (9) 

where  (-I)*”  accounts  for  the  change  in  direction  in  the  m^^  vertical 
leg,  and  the  Om  accounts  for  the  phase -change  from  one  leg  to 
another.  Let  us  say  that  the  distance  between  up  and  down  leg  is  N. 
Then  we  can  say 

E=Euexpi(ky-(ot)  +  Edexpi(-ky-o)t+C))  (10) 

where 

Eu  =  Eoexp-(z/ro)2,  Ed  =  Eoexp-((z-N)/ro)2  (11) 

In  Eqs.(10  and  11),  it  is  now  understood  that  E  is  in  the  x  direction,  z, 
the  distance  along  the  magnetic  field  on  which  the  beam  propagates 
now  replaces  the  transverse  variable  r. 

Expressing  the  electric  field  on  the  two  vertical  legs  in  the  form 
of  Eq.  (3  and  4),  and  Fourrier  decomposing  the  periodic  function  P, 
we  find 


E(y,z)  =  LnVitEo(ro/L)exp[-(k-i-qn)2ro2/4]expi(k-hqn)z 

(12) 

x{exp  iLky/A  +  exp  i[<D-(k+qn)P]exp-iLky/A}expi(kz-tot) 

where  qn  =  2%x\[L.  We  now  use  this  electric  field  in  the  linearized 
Vlasov  Equation  to  calculate  the  perturbed  distribution  function  and 
from  it,  the  perturbed  current.  The  calculation  reported  here  closely 
parallels  the  analysis  of  Vomvoridis  et  al^.  There  the  independent 
variables  were  the  guiding  center  variables,  related  to  the  normal 
Cartesian  variables  by 
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pjj=pj^cos(Qmz/pz  +<|)) 

Py  =plsin(Qmz/pz +<|)) 

X  =  Xg  +  (pj./nii2)sin(Qinz/pz +<t>)  (^) 

y  =  yg  .(pj_/mQ)cos(Qmz/pz  +<!>)  (d) 


(13) 


Then  the  linearized  Vlasov  Equation,  in  these  variables  becomes 
(9/3z  -i7ma)/pz)f  =  G(z) 


where 


G(z)  =  (YmeE(z)/(opz){2a)picos(Qmz7pz  +<|))9fo/3pi^} 


(15) 


where  fo  is  the  unperturbed  distribution  function  and  yis  the 
relativistic  factor  and  Q.  is  the  cyclotron  frequency.  In  Eq.(15)  we 
have  included  only  the  velocity  gradient  terms,  and  neglected  the 
terms  arising  from  gradients  in  Xg  and  yg.  These  typically  give  small 
corrections,  and  have  been  discussed  in  Ref.(lO).  The  perturbed 
distribution  function  f  is  then  given  by 


f  =  exp[iymtoz/pz]  Jzdz’exp[-iym(oz7pz]G(z') 


(16) 


From  the  perturbed  distribution  function,  we  calculate  the 
perturbed  current 


Jx  =-eJd3p  (px/mY)f 


(17) 


The  quantity  to  calculate  then  is  the  integral  of  Pi*Jx  over  the 
transverse  cross  section  and  over  a  periodicity  distance  L.  Without 
going  into  the  details  of  the  calculation  (it  is  quite  analogous  to  that 
in  Ref.(lO)),  we  find  that 

JPl*Jx  =  -Ajd^igd^pEn.s  [exp-(k+qn)7ro2/2]{2J’s2(Lkpi./ymQA) 


x2a)pi2/p2}9fo/api2/i[k+qn  +  (sftm  -  Ymto)/pz]  (18) 

where  A  =  Now  for  the  distribution  function,  take 

fo  =  I  {g(rg)mY/27cepz}  5(pi^“Pio^)  S(pz-Pzo) 
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where  I  is  the  beam  current  and  g(rg)  is  the  distribution  of  guiding 
centers,  normalized  so  that  |d2rgg(rg)=l.  Notice  that  the  distribution 
function  is  assumed  to  be  unperturbed  by  its  interaction  with  the 
radiation  upstream  from  the  periodicity  length  considered.  This  is 
examined  more  fully  in  the  appendix  where  it  is  shown  that  this 
interaction  (a  klystron  type  interaction)  has  a  negligible  effect,  at 
least  in  the  linear  regime. 

Consistent  with  our  assumptions  of  neglecting  the  various 
derivatives  with  respect  to  rg  in  the  derivation  of  Jx,  we  assume  that 
the  variation  of  g(rg)  is  small  over  distances  compared  to  k'l.  To  do 
the  integral  on  the  right  hand  side  of  Eq.(18),  we  perform  a  partial 
integration  over  pi^,  and  keep  only  those  terms  which  arise  from  a 
differential  of  the  resonant  denominator.  These  terms  are  the  ones 
which  have  that  denominator  raised  to  the  second  power,  and  are  the 
largest  terms  if  the  wave  frequency  is  nearly  equal  to  the  cyclotron 
harmonic.  To  continue,  we  approximate  the  effect  of  thermal  spread 
in  the  beam.  For  the  distribution  function,  let  us  assume  that  the 
beam  energy  (characterized  by  yor  the  magnitude  of  the  momentum 
p)  is  the  same  for  all  particles,  but  there  is  a  distribution  of  pitch 
angle  a  (pz/p  =  cos  a).  This  pitch  angle  spread,  assumed  to  be  small, 
is  characterized  by  5(x.  ^Ve  consider  the  effect  of  pitch  angle  spread 
only  in  the  resonant  denominator.  For  small  8a  this  is  the  only  place 
where  its  effect  can  be  significant.  The  integral  over  pz  can  be  done 
analytically  if  the  distribution  of  pitch  angle  is  approximated  by  a 
Lorentzian 


f(a)  =  {7r[(a-ao)2  +  (5a)2]}-l  (20) 

In  this  case  the  distribution  of  pz  is  also  Lorentzian  centered  at  pzo 
and  having  thermal  spread  pio5a.  While  we  cannot  use  the 
Lorentzian  to  accurately  calculate  the  effect  of  a  pitch  angle  spread,  it 
should  give  a  reasonable  extimate  of  when  the  spread  has  a 
significant  impact  on  the  growth  rate. 

If  we  assume  that  the  cyclotron  harmonics  are  all  decoupled, 
then  the  dispersion  relation  can  be  written  in  the  form 
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(21) 


-coA/cL  +k  =  Hs  [l-ipi8a/pz]-2  exp-0.5[(qn+k)ro]2 

x[k  +  qn  +  m(sa-Yco)/(pz-ipi5a)]-2 

where  to  simplify  the  notation,  in  Eq.  (21)  we  have  deleted  the  zero 
subscripts  on  the  momentum  variables.  Here, 

Hs  =  16Yecopi2l  J’s2(Lkpi/YmQA)  /(pz3  c3l2)  (22) 

This  then  summarizes  our  result  for  the  configuration  of  Fig  1.  For 
that  in  Fig  2,  an  analogous  calculation  gives  the  result 

-0)A/cL  +k  =  En  [Hs7sin2p]  [l-ipi6a/pz]-2  (23) 

xexp-0.5[(qn+k)ro/sin2p]2x[k  +Q)cosp/c  +  qn  +  m(sQ-YCo)/(pz-ipi5a)]-2 

where  p  is  the  angle  of  the  ray  path  to  the  horizontal  in  Fig  2  and 

Hs’  =  16Yeo)pi2l  J's2(Lkpisinp/YmaA)  / (pz3  c3l2)  (24) 

The  summation  in  |Pi*Jxis  a  summation  over  various  spatial 
harmonics  of  the  mode  structure,  each  spatial  harmonic  having  a 
different  resonant  denominator.  As  long  as  the  spatial  growth  rate 
Imk  is  much  less  than  2jt/L,  the  spatial  harmonics  can  be  considered 
one  at  a  time,  the  spatial  harmonic  with  k=-qn  being  the  dominant 
one.  Ultimately,  these  spatial  harmonics  become  small  due  to  the 
cutoff  of  the  Gaussian  exp-(k+qn)2ro2/2.  Hence  the  amplifier  can 
typically  operate  at  some  small  number  of  spatial  harmonics  before 
the  interaction  strength  becomes  negligibly  small.  Thus  the  fact  that 
the  system  is  periodic  with  periodicity  length  L  is  extremely 
important  in  that  it  breaks  the  eigenfunction  into  this  summation 
over  space  harmonics.  These  space  harmonics  are  sufficiently 
separated  in  wave  number  that  the  interaction  effectively  proceeds 
with  only  one  of  them  at  a  time.  This  greatly  increases  the  strength 
of  the  interaction  in  the  same  way  that  a  monoenergetic  beam 
interacts  more  strongly  with  a  wave  than  does  a  thermal 
distribution.  It  is  not  difficult  to  show  that  for  a  single  resonator,  the 
summation  over  k  becomes  an  integral,  and  the  conventional 
expression  (with  no  resonant  denominator)  is  obtained. 
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In  each  of  these,  the  assumption  is  that  the  radiation  mode  and 
beam  each  propagate  from  left  to  right.  This  is  the  usual  forward 
wave  amplifier  situation.  For  the  case  of  the  radiation  mode 
propagating  from  right  to  left,  the  backward  wave  case,  the 
dispersion  relations  are  the  same  except  that  the  sign  of  the  k  term 
on  the  left  hand  side  becomes  negative.  These  simple  dispersion 
relations  then  summarize  the  basic  theoretical  results  obtained  up  to 

now. 
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4.  Instability  condition 

In  a  traveling  wave  amplifier,  the  forward  wave  interaction 
cannot  be  absolutely  unstable,  but  the  backward  wave  can  be  We 
calculate  the  stability  threshold,,  paralleling  an  earlier  calculation! 2, 
We  only  consider  the  worst  case  of  zero  thermal  spread,  and  only  the 
fundamental  cyclotron  frequency.  Then  the  backward  wave 
dispersion  relations,  taking  only  a  single  spatial  harmonic,  derived  in 
the  last  section  all  had  the  form 

k(k-A)2  =  -Ko3  (25) 

where  A  and  Kq  are  both  real.  For  the  beam  energies  we  consider,  Kq 
decreases  as  a  function  of  cyclotron  harmonic  number.  As  we  will 
see,  the  stable  length  of  an  amplifier  is  expressed  as  coefficient  to  be 
derived,  divided  by  Kq.  Therefore,  all  cyclotron  harmonics  above  n 
will  be  stabilized  if  harmonic  n  is  itself  stabilized.  This  important 
result  for  the  quasi-optical  gyro  TWA  follows  directly  from  the  very 
simple  nature  of  the  spectrum,  especially  the  absence  of  cutoff 
modes.  In  other  gyro  TWA  experiments,  important  competing  modes 
have  turned  out  to  be  cutoff  modes  at  higher  harmonics. 

Consider  the  amplification  as  arising  from  the  coupling  of  a 
beam  mode  with  amplitude  Wi  and  a  radiation  mode  with  amplitude 
W2.  The  system  is  described  by  the  two  equations 


(k-A)2Wi  =  -Ko2W2 

(a) 

kW2  =  KoWi 

(b) 

This  is  cubic  system,  so  there  are  three  generally  complex  roots 
denoted  ki,k2and  ks.  The  coefficients  of  the  three  exponentials  are 
denoted  ^i,  ^2,  and  ^3.  Then  consider  a  system  of  length  'P.  The  beam 
mode  convects  from  left  to  right,  entering  at  z=0.  The  radiation  mode 
convects  from  right  to  left,  entering  the  system  at  z=T'.  The  condition 
for  instability  is  then  that  if  no  beam  mode  enters  at  0,  and  no 
radiation  mode  enters  at  'P,  there  will  be  a  radiation  mode  of 
amplitude  W20  exiting  the  system  at  z=0.  Since  the  beam  mode  is 
described  by  a  second  order  equation,  the  condition  for  no  beam 
mode  at  the  entrance  means  that  both  Wi  and  its  derivative  vanish 
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at  z=0.  Regarding  the  exponential  as  solutions  for  Wi,  we  find  that 
the  equations  for  the  three  mode  amplitudes  are 

0  =  iki^i+  ik2^2  +  ik3^3  (b) 

W20  =  (Ko/ki)^l+  (Ko/k2)^2  +  (Ko/k3)^3  (c) 

The  condition  for  instability  is  then  that  the  amplitude  of  W2  at  z='P 
vanishes,  or 

(^l/ki)expiki'P  +(^2/k2)expik2^  +(^3/k3)expik3'F  =  0  (28) 

Solving  for  the  coefficients,  the  instability  condition  becomes 

[(k3-k2)/ki]expiki'F  +  [(ki-k3)/k2]expik2'P 

+  [(k2-ki)/k3]expik3'P  =  0  (29) 

If  'F  is  very  small  compared  to  the  reciprocal  of  the  k’s,  which  in  turn 
scale  with  Kq,  clearly  there  is  no  solution  to  Eq.(29)  and  the  system  is 
stable.  As  'F  increases,  the  oscillatory  nature  of  the  exponentials  will 
allow  the  summation  to  equal  zero.  At  some  length  T'o.  there  will 
first  be  a  value  of  A=Ao  for  which  the  summation  vanishes.  At  this 
point  there  can  be  a  solution  of  a  non  zero  outgoing  radiation  mode 
with  no  inputs  of  either  radiation  or  beam  mode.  The  A©  then 
determines  the  frequency  of  the  instability,  and  the  determines 
the  maximum  length  (or  beam  current)  of  the  stable  system. 

A  numerical  solution  of  the  dispersion  relation  gives  the  result 


Ao=  -3.003/Ko 
'Fo  =  7.68/Ko 


(a) 

(b) 


(30) 
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5.  Results  and  Discussion 

Let  us  begin  by  considering  a  possible  configuration  for  both 
the  quasi-optical  and  corrugated  waveguide  configuration.  For 
normal  incidence,  if  the  mirror  -has  radius  of  curvature  Rm?  and  half 
the  separation  between  the  mirrors  is  Ly,  then  the  radiation  spot  size 
at  the  center  is  given  by  13 

ro  =  [?^Ly(Rm/Ly  -  l)l/2/7c]l/2  (31) 

where  X  is  the  radiation  wavelength.  The  spot  size  at  the  mirror  then 
turns  out  to  be 

tom  =  [XLy(Rm/Ly  -  l)-l/2/it]l/2  (32) 

If  we  take  a  spot  size  to  of  0.5cm,  a  wavelength  of  0.3  cm,  and  Ly  of  3 
cm,  we  find  that  the  mirror  radius  Rm  is  about  5  cm  and  the  spot  size 
at  the  mirror  is  about  .62  cm.  To  reduce  diffraction  loss,  the  mirror 
size  should  be  a  radius  of  about  1.5  cm.  Since  the  mirror  is  tilted  at 
45°,  the  total  half  length  of  the  mirror  in  the  vertical  direction  is 
about  2  cm,  so  the  vertical  ray  path  is  about  10  cm.  Hence  a  magnet 
bore  of  roughly  4  inches  should  be  sufficient  to  hold  the  tube.  If  the 
facing  mirror  edges  are  two  centimeter  apart,  the  horizontal  length  of 
the  ray  path  is  L=10  cm,  and  the  total  ray  path  within  the  periodicity 

cell  is  A  =  22  cm. 

The  corrugated  waveguide  propagates  TEM  modes  with  Ex 
proportional  to  Jo(2.4r/a),  with  very  low  loss.  If  this  mode  is 
decomposed  into  a  summation  of  Gaussian  modes  with  ro  =  0.64a, 
then  about  98%  of  the  poower  is  in  the  fundamental  Gaussian 
models.  Thus  the  Gaussian  mode  is  an  excellent  approximation  to 
the  field  in  the  corrugated  waveguide  as  well.  For  this  case  an 
important  consideration  is  the  loss  on  transiting  the  gap  and  on 
reflection  at  the  mitre  bends.  According  to  Doane  and  Moeller  5, 
the  loss  on  transiting  the  gap  is  given  by 


Gap  Loss  (db)  =  1.7  [QX/2a2]3/2 
Mitre  Bend  Loss  (db)  =  0.85  [>./a]3/2 


(33) 
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where  a  is  the  waveguide  radius  and  Q  is  the  gap  space.  For  our  case 
of  X  =0.3  cm,  a  waveguide  radius  of  1 .5  cm  and  a  gap  length  of  also 
1.5  cm,  we  find  that  there  is  a  1.2%  loss  at  each  gap  and  a  1.7%  loss  at 
each  mitre  bend.  Since  there  are  2  gaps  and  4  mitre  bends  in  each 
periodicity  length  L,  there  is  about  a  9%  loss  in  the  periodicity  cell. 
This  loss  is  only  loss  from  the  TEM  mode;  much  of  the  lost  power  is 
still  bouncing  alfound  the  waveguide.  If  the  wave  amplifies  strongly 
from  one  periodicity  cell  to  the  next,  the  output  radiation  should  be 
about  90%  in  the  desired  mode. 

Now  let  us  evaluate  the  growth  and  bandwidth  from  the 
dispersion  relations  calculated  in  Section _ 3.  We  consider  only  the 
case  where  the  space  harmonics  decouple  from  one  another,  so  that 
the  beam  interacts  with  only  one  of  them.  The  use  of  a  single  spatial 
harmonic  itself  puts  a  limit  on  the  bandwidth  which  can  be 
accurately  calculated.  Since  k  =  coA/cL,  and  a  change  in  k  of  27t/L 
brings  us  to  the  next  spatial  harmonic,  the  maximum  accurate 
bandwidth  we  can  calculate  with  a  single  spatial  harmonic  is  about 
1.5  GHz  for  the  parameters  specified.  1.5  GHz. 

If  we  define 


K  =  k  -  coA/cL 

and  assume,  strictly  for  simplicity  sake  that  the  cyclotron  frequency 
corresponds  exactly  to  cLqn/A,  and  define  co  =  sf2/Y+5co,  where  s  is  the 
cyclotron  harmonic,  we  find  that  the  dispersion  relation  for  a  cold 
beam  is 


K[K-A]2=Ko3  (34) 

If  both  K  and  A  are  scaled  with  Kq,  then  K/Kq  is  a  function  of  only  the 
single  parameter  A/Kq.  A  numerical  solution  for  the  fastest  growing 
root  is  shown  in  Fig. (8).  The  minimum  growth  length  turns  out  to  be 
about  1.2/Ko. 

Let  us  consider  an  interaction  at  the  cyclotron  frequency  for  a 
configuration  like  that  shown  in  Fig.(l).  In  this  case,  we  find  that 
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Ko^  =  Kylepi^coJ/Lpz^c^L^] 

A  =  [A/cL  -Ym/(pz-ipi5a)]8co 


(a) 


(35) 

(b) 

where  we  have  actually  included  the  effect  of  the  thermal  spread  in 
the  resonant  denominator.  Let  us  consider  a  particular  beam  for 
this  configuration.  Take  an  energy  of  100  kV,  a  value  of  tan  a  of  1.5 
and  no  thermal  spread.  This  is  a  fairly  high  value  for  an  amplifier 
configuration.  However  our  hope  is  that  this  configuration  will  be 
much  more  stable  than  a  waveguide  configuration,  for  which  such  a 
high  value  of  pi/pz  would  almost  surely  be  unstable.  Then  if  © 
=6x10^^  s'l,  we  find  that 

Ko3  (cm-3)  =  3.2x10-2  [l(Amps)/L2(cm2)]  (36) 

Thus  if  we  consider  a  beam  of  ten  amps  and  L=10  cm,  the  minimum 
growth  length  is  about  8  cm,  and  maximum  stable  length  is  about  60 
cm.  One  real  advantage  of  the  configuration  of  Fig.(l)  is  that  it  is  not 
very  sensitive  to  beam  thermal  spread.  As  is  apparent  from 
Eq.(35b),  the  condition  for  neglect  of  thermal  spread  at  the  band 

center  (5©  =  0)  is 

Pz  »  pi6a  (36) 

a  condition  not  very  difficult  to  satisfy  for  typical  gyrotron  beams. 

As  5©  increases  however,  the  effect  of  thermal  spread  becomes  more 
pronounced  because  of  its  contribution  in  the  resonant  denominator. 
As  a  function  of  5©  the  condition  for  the  neglect  of  thermal  spread  is 

5a  5©/vz  «  Ko  (37) 

For  the  beam  parameters  we  have  chosen,  and  5©  chosen  to  be  2k 
times  1.5  GHz,  we  see  that  the  maximum  allowed  angular  spread  is 
about  10*1  radians. 

Let  us  now  calculate  the  bandwidth,  for  the  cold  beam,  in  the 
linear  regime.  Imagine  that  the  input  power  is  100  W  and  the 
output  power  is  10^  W,  so  the  power  gain  is  a  thousand.  This  is 
about  3.5  e  folds  in  amplitude.  At  the  band  edge,  the  power  is  down 
by  a  factor  of  2,  or  the  amplitude  by  factor  of  1.4,  (ie  expO.3).  From 
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the  graph  in  Fig. (8),  this  gives  a  spread  in  A/Ko  of  about  2.  For  the 
parameters  chosen  this  means  a  bandwidth  of  about  1.5  GHz,  the 
maximum  accurate  bandwidth.  Thus  a  more  accurate  calculation  of 
the  bandwidth  must  involve  the  calculation  of  the  simultaneous 
interaction  of  several  spatial  harmonics  and  thermal  spread. 

At  the  cyclotron  harmonic,  for  the  same  system  except  now  a 
beam  energy  of  150  kV,  we  find  that  Kq^  is  smaller  by  about  a  factor 
of  0.28,  so  the  minimum  growth  length  now  turns  out  to  be  about  12 
cm.  The  total  length  must  still  be  below  that  determined  for  the 
cyclotron  frequency  instability  however,  or  else  severs  or  directional 
couplers  would  have  to  be  used.  . 

Finally,  let  us  consider  the  configuration  of  Fig.(2).  Here, 
interaction  strengths  and  bandwidths  are  somewhat  larger,  and 
required  magnetic  fields  are  somewhat  lower,  because  of  the  Doppler 
shift  provided  by  the  beam's  velocity  parallel  to  the  direction  of 
propagation.  However  this  seems  to  be  counterbalanced  by  the  fact 
that  the  beam  requirements  are  much  more  difficult  to  achieve. 

From  Eq.(23),  the  requirements  on  beam  thermal  spread  are  given 
by 

Kpz  » pi5a(a)/c)cosp  (38) 

This  is  a  much  more  difficult  condition  to  satisfy.  For  the  beam 
parameters  we  have  been  taking,  this  would  mean  an  a  spread  of 
about  a  percent  or  less. 
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Appendix  The  Effect  of  Prebunching  in  the  Linear  Regime 

Here,  we  consider  the  effect  on  Pi*Jxfrom  the  the  upstream 
periodicity  elements.  That  is  in  Sec.  3,  when  calculating  Pj  Jx,  we 
assumed  that  the  input  distribution  function  was  a  delta  function 
P±,  Pz,  and  had  uniform  distribution  in  azimuthal  angle  4).  However 
actually,  the  beam,  by  the  time  it  enters  the  periodicity  unit  in 
ruction,  has  passed  through  all  of  the  periodicity  units  upstream. 
Thus  the  distribution  function  is  actually  not  a  delta  function,  but 
must  also  reflect  the  interaction  with  the  upstream  cavities.  ^ 

consider  the  effect  of  this  interaction  in  the  linear  regime.  We  show 
that  the  effect  is  very  small  compared  to  those  retained  in  Sec.3. 

We  only  consider  the  effect  from  the  nearest  upstream 
periodicity  unit,  and  also  consider  only  the  effect  of  prebunching 
from  that  part  of  the  radiation  going  in  the  same  direction  in  y  m 
each  unit.  As  we  will  see,  there  are  phase  averages,  which  nearly 
vanish,  involved  in  the  integral  of  Pi  *Jx  calculated  on  the  radmtion 
going  up  in  y,  resulting  from  the  prebunching  of  the  electron  beam, 
in  the  nearest  upstream  periodicity  unit,  from  the  radiation  going 
down  in  y. 

One  can  integrate  the  linearized  equation  for  pi  through  the 
upstream  radiation  field  and  find  that  the  wave  going  in  the  positive 
direction  in  y  gives  the  the  i*  electron  of  the  beam  a  change  in  pi 

given  by 

Apii=  N7teEoro/2mvz]  cos(a)yi/c  -  coti  +<^i)  exp-[Q/Y  -0)]2ro2/(2vz)2  (Al) 

where  y\  is  the  guiding  center  position  of  the  entering  electron,  q  is 
the  entrance  time,  and  <l)i  is  the  entrance  gyro  phase. 

The  change  in  Apii  gives  a  change  to  the  y  of  each  electron, 
which  in  turn  gives  a  change  to  its  orbital  frequency  between  the  n-1 
and  n  periodicity  units.  When  it  enters  the  next  periodicity  unit,  the 

phase  angle  of  the  i^^  electron  is 


(1)  =  4)1  +  QL/yvz  -  OLpiApii/2Y3vzm2c^ 


(A2) 
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Then,  in  the  integral  over  t  and  ^  integral  for  Pi*Jx,  we  transform 
these  variables  to  t’  and  0’,  and  find  that  the  linear  expresion  is 
given  by 

/d3rd3pdt’  Pi*Jx  =  Jd3rd3pdt’We6s(ky+<t>i-coti)sin((t)i+nz/YVz) 

xcos[ky-a)(ti+z/vz)  +^]  exp-[(z+L)/ro]2  (A3) 

where  C  is  the  phase  shift  in  going  across  the  periodicity  unit  and 


W  =  V7c[QLe2EroPi2fo(pi2,pz,(|)i)/4m3vz273] 

X  exp-[(Q/Y  -  a)]2(ro/2vz)2 


(A4) 


Here,  fo  is  assumed  to  be  distributed  uniformly  in  <t)i.  The 
trigonometric  integral  in  Eq.(A4)  breaks  up  into  several  terms  which 
integrate  to  zero  over  phase.  The  one  that  does  not  is 

Jd3rd3pdt’cos2(ky+<t)i-a)ti)  sin^  sin[(Q/Y-co)z/vz] 

Note  that  had  the  downward  traveling  wave  provided  the 
prebunching  and  the  upward  wave  the  extraction,  the  first  factor 
above  would  have  been,  cos(ky+(}>i-(oti)cos(-ky+(l)i-(oti),  which  would 
also  integrate  to  zero  over  phase.  Doing  the  phase  integrals  Eq.(A3), 
we  find  that 


Jd3rd3pdt’  Pi*Jx=W  sinC  sin[(a/Y-co)L/vz] 
X  exp-2[(D/Y  -  0)]2(ro/2vz)^ 


(A5) 


Thus  in  this  contribution  to  the  energy  input,  there  is  no  resonant 
denominator  rasied  any  power.  Hence  this  term  is  but  a  small 
correction  to  those  calculated  in  Sec.(3)  which  typically  had 
denominators  like  [(D/y  -  co]‘2. 
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Schematic  of  the  oblique  incidence  quasi-optical  gyro  TWA. 
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CATHODE^ 

Schematic  of  the  quasi-optical  gyrotron  and  gyro  klystron. 


Schematic  of  the  Lincoln  Laboratory  clam  shell  reflector. 


Schematic  of  the  Lincoln  Laboratory  quasi-optical  circulator. 
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8.  Solution  of  the  dispersion  relation  for  K  as  a  function  of  A. 
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